A measure of the non- Gaussian character of a quantum state 
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We address the issue of quantifying the non-Gaussian character of a bosonic quantum state and introduce 
a non-Gaussianity measure based on the Hilbert-Schmidt distance between the state under examination and a 
reference Gaussian state. We analyze in details the properties of the proposed measure and exploit it to evaluate 
the non-Gaussianity of some relevant single- and multi-mode quantum states. The evolution of non-Gaussianity 
is also analyzed for quantum states undergoing the processes of Gaussification by loss and de-Gaussification by 
photon-subtraction. The suggested measure is easily computable for any state of a bosonic system and allows 
to define a corresponding measure for the non-Gaussian character of a quantum operation. 

PACS numbers: 03.67.-a, 03.65.Bz, 42.50.Dv 



I. INTRODUCTION 

Gaussian states play a crucial role in quantum information 
processing with continuous variables. This is especially true 
for quantum optical implementations since radiation at ther- 
mal equilibrium, including the vacuum state, is itself a Gaus- 
sian state and most of the Hamiltonians achievable within the 
current technology are at most bilinear in the field operators, 
i.e. preserve the Gaussian character |H 0. Ht]. As a matter 
of fact, using single-mode and entangled Gaussian states, lin- 
ear optical circuits and Gaussian operations, like homodyne 
detection, several quantum information protocols have been 
implemented, including teleportation, dense coding and quan- 
tum cloning |01 . 

On the other hand quantum information protocols required 
for long distance communication, as for example entangle- 
ment distillation and entanglement swapping, rely on non- 
Gaussian operations. In addition, it has been demonstrated 
that teleportation JH 0, 01 and cloning |H] of quantum states 
may be improved by using non-Gaussian states and non- 
Gaussian operations. Indeed, de-Gaussification protocols for 
single-mode and two-mode states have been proposed 1111^,01 
and realized jgt]. It should be also noticed that any strongly 
superadditive function is minimized, at fixed covariance ma- 
trix, by Gaussian states. This is crucial to prove extremality 
of Gaussian states and Gaussian operations ifiol [Till f or what 
concerns various quantities as channel capacities lll2ll . multi- 
partite entanglement measures 11311 and distillable secret key 
in quantum key distribution protocols. Since in most cases 
these quantities can be computed only for Gaussian states, a 
non-Gaussianity measure may serve as a guideline to quan- 
tify them for the class of non-Gaussian states. Overall, non- 
Gaussianity is revealing itself as a resource for continuous 
variable quantum information, and thus we urge a measure 
able to quantify the non-Gaussian character of a quantum 
state. 

In this paper we introduce a novel quantity, the non- 



Gaussianity 8[q] of a quantum state, which quantifies how 
much a state fails to be Gaussian. Our measure, which is based 
on the Hilbert-Schmidt distance between the state itself and a 
reference Gaussian state, can be easily computed for any state, 
either single-mode or multi-mode. 

The paper is structured as follows. In the next Section we 
introduce notation and review the basic properties of Gaussian 
states. Then, in Section[nI]we introduce the formal definition 
of 8[q] and study its properties in details. In Section [Tvl we 
evaluate non-Gaussianity of relevant quantum states whereas 
in Section |V1 we analyze the evolution of non-Gaussianity for 
known Gaussification and de-Gaussification maps. SectionlVll 
closes the paper with some concluding remarks. 



H. GAUSSIAN STATES 

For concreteness, we will use here the quantum optical ter- 
minology of modes carrying photons, but our theory applies 
to general bosonic systems. Let us consider a system of n 
modes described by mode operators Ofc, k = 1 . . . n, satis- 
fying the commutation relations [cifc,aj] = 8k j. A quantum 
state q of the n modes is fully described by its characteristic 
function H 



X [q](\)=Tt[qD(\)] 

where D(X) = (£fZ=i Dk(Xk) is the n-mode displacement 
operator, with A = (Ai, . . . , A„) T , £ (D, and where 

D k {\k) = cxp{X k a f k - X* k a k } 

is the single-mode displacement operator. The canonical op- 
erators are given by: 
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with commutation relations given by [qj,p k ] = iSjk- Upon in- 
troducing the real vector R = (qx,pi, . . . , q n ,Pn) T , the com- 
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mutation relations rewrite as 

[R k ,R j }=ifl kj 

where Q;y are the elements of the symplectic matrix ft = 
i ©)! =1 (7 2? <J2 being the y-Pauli matrix. The covariance ma- 
trix <r = cr[g] and the vector of mean values X = X[g] of a 
quantum state g are defined as 

Xj = (Rj) 

v kj = ±({R k ,R j })-(R j )(R k ) 

(1) 

where {A, B} = AB + BA denotes the anti-commutator, and 
(O) = Tr[g O] is the expectation value of the operator O. 

A quantum state go is referred to as a Gaussian state if its 
characteristic function has the Gaussian form 

x[qg](A) = exp |-^A T tA + X T OA| 

where A is the real vector A = 
(ReAi,ImAi, . . . ,ReA„,ImA n ) T . Of course, once the 
covariance matrix and the vector of mean values are given, a 
Gaussian state is fully determined. For a single-mode system 
the most general Gaussian state can be written as 

Q G = D(a)S(CMn t )&(QDH<x), 

D(a) being the displacement operator, S(C) = 
exp[iC(a t ) 2 — \C,* a2 \ m e squeezing operator, a,(sC, and 
v(n t ) = (1 + n t ) _1 [n t /(l + n t )] a ' a a thermal state with n t 
average number of photons. 

III. A MEASURE OF THE NON-GAUSSIAN CHARACTER 
OF A QUANTUM STATE 

In order to quantify the non-Gaussian character of a quan- 
tum state g we use a quantity based on the distance between 
g and a reference Gaussian state r, which itself depends on g. 
Specifically, we define the non-Gaussianity 5[g] of the state g 
as 

s\e] - (2) 

where Dhs[q,t] denotes the Hilbert-Schmidt distance be- 
tween g and r 

r>2 r 1 ^t/ n2i M[gj +MM ~ 2«[g,r] 

= 2 TT ^ g - T ) 1 = 2 ' ( ) 

with = Tr[g 2 ] and r] = Tr[gr] denoting the purity of 
g and the overlap between g and r respectively. The Gaussian 
reference r is the Gaussian state such that 

X[g] = X[t] 
tr[g] = <t[t] 



i.e. t is the Gaussian state with the same covariance matrix a 
and the same vector X of the state g. 

The relevant properties of 6[g], which confirm that it repre- 
sents a good measure of the non-Gaussian character of g, are 
summarized by the following Lemmas: 
Lemma 1: S[g] — iff g is a Gaussian state. 
Proof: If S[g] = then g = r and thus it is a Gaussian state. 
If g is a Gaussian state, then it is uniquely identified by its first 
and second moments and thus the reference Gaussian state r 
is given by r = g, which, in turn, leads to Dhs [q, t] = and 
thus to 6[g] = 0. 

Lemma 2: If U is a unitary map corresponding to a symplec- 
tic transformation in the phase space, i.e. if U = exp{— iH} 
with hermitian i7 that is at most bilinear in the field operators, 
then 6[UqW] = S[g]. This property ensures that displace- 
ment and squeezing operations do not change the Gaussian 
character of a quantum state. 

Proof: Let us consider g' = UgW . Then the covariance ma- 
trix transforms as cr[g'} = Eer[g]£ T , £ being the symplectic 
transformation associated to U. At the same time the vector of 
mean values simply translates to X' = X + Xq, where Xq 
is the displacement generated by U. Since any Gaussian state 
is fully characterized by its first and second moments, then the 
reference state must necessarily transform as r' = UtU', i.e. 
with the same unitary transformation U. Since the Hilbert- 
Schmidt distance and the purity of a quantum state are invari- 
ant under unitary transformations the lemma is proved. 
Lemma 3: 5[g] is proportional to the squared L 2 (<C n ) dis- 
tance between the characteristic functions of g and of the ref- 
erence Gaussian state r. In formula: 

S[g\ oc J d 2n A[xb](A)-xH(A)] 2 . (4) 

Since the notion of Gaussianity of a quantum state is de- 
fined through the shape of its characteristic function, and since 
the characteristic function of a quantum state belongs to the 
L 2 ((D") space Hill , we address L 2 (<C) distance to as a good 
indicator for the non Gaussian character of g. 
Proof: Since characteristic functions of self-adjoint operators 
are even functions of A and by means of the identity 

/d 2n \ 
±— x [0 1 }(X) X [0 2 }(-\), 

we obtain 

Lemma 4: Consider a bipartite state g = g& ® g G . If go is a 
Gaussian state then S[g] = 6[qa\- 
Proof: we have 

iAq\ = iA.qa]iAqg\ 

h[t] = h[t a ][i{t g } 
k[q, t] = k[qa,ta]k[qg, Qg] ■ 
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Therefore, since k[qc, Qg\ = iA.Qg\ we arrive at 

V>[qaMqg] + v[ta]v[qg] ~ 2k[q a , t a ]k[q g , q g ] 



S[q] 



2^i[g A ]^[g c 



(5) 



The four properties illustrated by the above lemmas are the 
natural properties required for a good measure of the non- 
Gaussian character of a quantum state. Notice that by using 
the trace distance [g, t] = iTr| g— t\ instead of the Hilbert- 
Schmidt distance we would lose Lemmas 3 and 4, and that the 
invariance expressed by Lemma 4 holds thanks to the renor- 
malization of the Hilbert-Schmidt distance through the purity 
fi[g\. We stress the fact that our measure of non-Gaussianity 
is a computable one: It may be evaluated for any quantum 
state of n modes by the calculation of the first two moments 
of the state, followed by the evaluation of the overlap with the 
corresponding Gaussian state. 

Notice that S[g] is not additive (nor multiplicative) with re- 
spect to the tensor product. If we consider a (separable) multi- 
partite quantum state in the product form g = ® r k L =1 gk, the 
non-Gaussianity is given by 



S[Q] = 



2 nLiM£'fc] 



(6) 



where is the Gaussian state with the same moments of gk- 
In fact, since the state g is factorisable, we have that the cor- 
responding Gaussian r is a factorisable state too. 



IV. NON-GAUSSIANITY OF RELEVANT QUANTUM 
STATES 

Let us now exploit the definition © to evaluate the non- 
Gaussianity of some relevant quantum states. At first we con- 
sider Fock number states \p) of a single mode as well as mul- 
timode factorisable states \p)® n made of n copies of a num- 
ber state. The reference Gaussian states are a thermal state 
T. p = v(p) with average photon number p and a factorisable 
thermal state tn = [^(p)]®" with average photon number p 
in each mode lfl5ll . Non-Gaussianity may be analytically eval- 
uated, leading to 



%»<Pll = i(i + ^ 



1 



p+1 \P+1 



5[(\P)(P\) 



0)1] 



1 



1 



2p+l 



p+1 \p+l 



In the multimode case of |p}®", we seek for the number of 
copies that maximizes the non-Gaussianity. In Fig. [JJ we 
show both S p = 6[\p)(p\] and 5 p = max n S[(\p) (p|)® n ] as a 
function of p. As it is apparent from the plot non-Gaussianity 
of Fock states \p) increases monotonically with the number 
of photon p with the limiting value S p = 1/2 obtained for 
p — > oo. Upon considering multi-mode copies of Fock states 
we obtain larger value of non-Gaussianity: S p is a decreasing 
function of p, approaching S — 1/2 from above. The value 




FIG. 1: (Top): Non-Gaussianity of single mode Fock states (gray) 
\p) and of multi-mode Fock states |p)® n (black) as a function of p. 
Non-Gaussianity for multi-mode states has been maximized over the 
number of copies n. (Bottom): Non-Gaussianity, as a function of the 
parameter <j>, for the two-mode superpositions |$)) (dashed gray), 
|*1/)) (solid gray), and for the single-mode superposition of coherent 
states \ips) for a — 0.5 (solid black) and a — 5 (dashed black). 

of S p corresponds to n = 3 for p < 26 and to n — 2 for 
27 < p < 250. 

Another example is the superposition of coherent states 

H's) = AT' 1 /2 (cos c/)\a) +sin<£| -a)) (7) 

with normalization TV = 1 + sin(2</>) exp{— 2a 2 } which for 
4> = ±7r/4 reduces to the so-called Schrodinger cat states, 
and whose reference Gaussian state is a displaced squeezed 
thermal state r s = D(C)S(r)v(N)S^ (r)D^ (C), where the 
real parameters C, r, and N are analytical functions of <f> and 
a. Finally we evaluate the non-Gaussianity of the two-mode 
Bell-like superpositions of Fock states 



|$)) = cos</>|0,0) 
|*)) = cos<£|0,l) 



sin < 



.11,1) 
sin</>|l, 0), 



which for <ft = ±7r/4 reduces to the Bell states |$ ± ) 
and | V I' ± ). The corresponding reference Gaussian states 
are respectively a two mode squeezed thermal state r$ = 
S 2 (0HN) ® v(N)]Sl(0, where S 2 (Q = exp^a\4 - 
£,*ab) denotes the two-mode squeezing operator, and = 
R(9)[v(Ni)® v(N 2 )]R' ! (9), namely the correlated two-mode 
state obtained by mixing a single-mode thermal state with 
the vacuum at a beam splitter of transmissivity cos 2 6, i.e. 
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R(0) = exp[i9(a[a2 + a 2 ai)]. All the parameters involved in 
these reference Gaussian states are analytical functions of the 
superposition parameter (j>. Non-Gaussianities are thus evalu- 
ated by means of (fJI and are reported in Fig. Q]as a function 
of the parameter <fi. As it is apparent from the plot, the non- 
Gaussianity of single-mode states does not surpass the value 
S = 1/2, and this fact is confirmed by other examples not 
reported here. 

As concern the cat-like states, we notice that for small val- 
ues of a the non-Gaussianity of the superposition \ips) shows 
a different behavior for positive and negative values of the pa- 
rameter (j>: for <j> > and a = 0.5 we have almost zero S, 
while higher values are achieved for <f> < 0. For higher val- 
ues of a (a = 5 in Fig. Q]), non-Gaussianity becomes an even 
function of <f>. This different behavior can be understood by 
looking at the Wigner functions of even and odd Schrodinger 
cat states for different values of a: for small values of a the 
even cat's Wigner function is similar to a Gaussian function, 
while the odd cat's Wigner function shows a non-Gaussian 
hole in the origin of the phase space; increasing the value of a 
the Wigner functions of the two kind of states become similar 
and deviate from a Gaussian function. 

We have also done a numerical analysis of non-Gaussianity 
of single-mode quantum states represented by finite superpo- 
sition of Fock states 

d 

Qd = ^2 ^k\n){k\ . (8) 

n,k=0 

To this aim we generate randomly quantum states in a finite 
dimensional subspaces, dim(_ff ) = d + 1 < 21, following the 
algorithm proposed by Zyczkowski et al flflffft , i.e. by gen- 
erating a random diagonal state (i.e. a point on the simplex) 
and a random unitary matrix according to the Haar measure. 
In Fig. |2]we report the distribution of non-Gaussianity S[gd), 
as evaluated for 10 5 random quantum states, for three different 
value of the maximum number of photons d. As it is apparent 
from the plots the distribution of S[gd] becomes Gaussian-like 
for increasing d. In the fourth panel of Fig. [2] we thus re- 
port the mean values and variances of the the distributions as 
a function of the maximum number of photons d. The mean 
value increases with the dimension whereas the variance is a 
mono tonic ally decreasing function of d. 

Also for finite superpositions simulations we did not ob- 
serve non-Gaussianity higher than 1/2. Therefore, although 
we have no proof, we conjecture that 8 = 1/2 is a limiting 
value for the non-Gaussianity of a single-mode state. Higher 
values are achievable for two-mode or multi-mode quantum 
states (e.g. 5 = 2/3 for the Bell states I**))). 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 



<5[p 2 ] 5 [Pio] 




<5[P 20 ] d 

FIG. 2: Distribution of non-Gaussianity S[gd] as evaluated for 10 5 
random quantum states, for three different value of the maximum 
number of photons d. Top: d = 2 (left), d = 10 (right); Bottom: 
d = 20 (left). (Bottom-right): Mean values and variances of the non- 
Gaussianities evaluated for 10 J random quantum states, as a function 
of the maximum number of photons d. 



with a bath of oscillators at zero temperature. This is per- 
haps the simplest example of a Gaussification protocol. In 
fact the interaction drives asymptotically any quantum state to 
the vacuum state of the harmonic system, which, in turn, is a 
Gaussian state. The evolution of the system is governed by the 
Lindblad Master equation g = %£[a]g, where g denotes time 
derivative, 7 is the damping factor and the Lindblad superop- 
erator acts as follows C[a]g = 2al ga — a'ag — ga)a. Upon 
writing r} = e~ 7 * the solution of the Master equation can be 
written as 

g(v) = Y, v ™e v rn ( 9 ) 

m 

V m = [(l~ V ) m /m^a m r 1 ^ aia - m \ 

where g is the initial state. In particular for the system ini- 
tially prepared in a Fock state g p = \p)(p\, we obtain, after 
evolution, the mixed state 



V. GAUSSIFICATION AND DE-GAUSSIFICATION 
PROCESSES 

We have also studied the evolution of non-Gaussianity 
of quantum states undergoing either Gaussification or de- 
Gaussification processes. First we have considered the Gaus- 
sification of Fock states due do the interaction of the system 



v 

m ;=o 

witha/. p (?7) = (^) (1 — rj) p ^ l ij l . The reference Gaussian state 
corresponding to g p (r]) is a thermal state t p (ij) = v(pr)) with 
average photon number pi]. Non-Gaussianity of g p (r]) can be 
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evaluated analytically, we have 

Spr, = 6[q p {ti)] 



1 



2(1 - r?) 2 ™ 2 Fx (-m,-m,l;^ 
x ^(l- V ) 2m 2 F 1 (-m,-m,l; 

+ (1 + 2mr\) 



(V I) 2 
_! 2(1 + (m - I)??)' 



(1 + m7j) m+1 



(11) 



2-Fi(a, 6, c; x) being a hypergeometric function. We show the 
behavior of 5 prj in Fig. [3] as a function of 1 — 77 for different 
values of p. As it is apparent from the plot <5™ is a monotoni- 
cally decreasing function of 1 — 77 as well as a monotonically 
increasing function of p. That is, at fixed time t the higher 
is the initial photon number p, the larger is the resulting non- 
Gaussianity. 




0.6 0.8 




FIG. 3: (Left): Non-Gaussianity of Fock states \p) undergoing Gaus- 
sification by loss mechanism due to the interaction with a bath of os- 
cillators at zero temperature. We show S vp as a function of 1 — 77 
for different values of p: from bottom to top p = 1, 10, 100, 1000. 
(Right): Non-Gaussianity of gips as a function of T for r — 0.5 and 
for different values of e — 0.2, 0.4, 0.6, 0.8 (from bottom to top). 
Sips results to be a monotonous increasing function of T, while e 
only slightly changes the non-Gaussian character of the state. 



Let us now consider the de-Gaussification protocol ob- 
tained by the process of photon subtraction. Inconclusive Pho- 
ton Subtraction (IPS) has been introduced for single-mode and 
two-mode states in |@, 0, [jjl] and experimentally realized in 
J^]. In the IPS protocol an input state g( m ) is mixed with 
the vacuum at a beam splitter (BS) with transmissivity T and 
then, on/off photodetection with quantum efficiency e is per- 
formed on the reflected beam. The process can be thus charac- 
terized by two parameters: the transmissivity T and the detec- 
tor efficiency e. Since the detector can only discriminate the 
presence from the absence of light, this measurement is in- 
conclusive, namely it does not resolve the number of detected 
photons. When the detector clicks, an unknown number of 
photons is subtracted from the initial state and we obtain the 
conditional IPS state gips- The conditional map induced by 
the measurement is non-Gaussian 01 > an d me output state is 
de-Gaussified. Upon applying the IPS protocol to the (Gaus- 
sian) single-mode squeezed vacuum S(r)|0) (r £ R), where 
S(r) is the real squeezing operation we obtain lfl8ll the con- 
ditional state qips, whose characteristic function x[qips]W 
is a sum of two Gaussian functions and therefore is no longer 
Gaussian. The corresponding Gaussian reference state is a 
squeezed thermal state tips = S(^ips)v(Nips)S^ (£ips) 



where the parameters £/ps and Nips are analytic functions of 
r, T and e. Non-Gaussianity Sips = <5/ps(T, e, r) has been 
evaluated, and in Fig. [3] (right) we report Sips for r = 0.5 
as a function of the transmittivity T for different values of 
the quantum efficiency e. As it is apparent from the plot the 
IPS protocol indeed de-Gaussifies the input state, i.e. nonzero 
values of the non-Gaussianity are obtained. We found that 
Sips is an increasing function of the transmissivity T which 
is the relevant parameter, while the quantum efficiency e only 
slightly affects the non-Gaussian character of the output state. 
The highest value of non-Gaussianity is achieved in the limit 
of unit transmissivity and unit quantum efficiency 

Vim 6 IPS =8[\l){l\] = 6[S{r)\l)(l\&(r)], 

T,7?— >1 

where the last equality is derived from Lemma 2. This result 
is in agreement with the fact that a squeezed vacuum state 
undergoing the IPS protocol is driven towards the target state 
S(r)\l) in the limit of T, e — > 1 OH. Finally, we notice that 
for T, e 7^ 1 and for r — > 00 the non-Gaussianity vanishes. In 
turn, this corresponds to the fact that one of the coefficients 
of the two Gaussians of x[qips]W vanishes, i.e. the output 
state is again a Gaussian one. 



VI. CONCLUSION AND OUTLOOKS 

Having at disposal a good measure of non-Gaussianity for 
quantum state allows us to define a measure of the non- 
Gaussian character of a quantum operation. Let us denote 
by Q the whole set of Gaussian states. A convenient defi- 
nition for the non-Gaussianity of a map 8 reads as follows 
6[£] = max e6 g S[£(g)], where £(g) denotes the quantum 
state obtained after the evolution imposed by the map. Indeed, 
for a Gaussian map £ g , which transforms any input Gaussian 
state into a Gaussian state, we have S[£ g ] = 0. Work along 
this line is in progress and results will be reported elsewhere. 

In conclusion, we have proposed a measure of the non- 
Gaussian character of a CV quantum state. We have shown 
that our measure satisfies the natural properties expected from 
a good measure of non-Gaussianity, and have evaluated the 
non-Gaussianity of some relevant states, in particular of states 
undergoing Gaussification and de-Gaussification protocols. 
Using our measure an analogue non-Gaussianity measure for 
quantum operations may be introduced. 
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APPENDIX A: GAUSSIAN REFERENCE WITH 
UNCONSTRAINED MEAN VALUE 

As we have seen from the above examples S[g] of Eq. (f2j) 
represents a good measure of the non-Gaussian character of a 
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quantum state. A question arises on whether different choices 
for the reference Gaussian state r may lead to alternative, 
valid, definitions. As for example (for single-mode states) we 
may define 



minB 2 HS [g,r]/fx[g], 



(Al) 



where r = D(C)S(£)v(N)S* (QD* (C) is a Gaussian state 
with the same covariance matrix of g and unconstrained vec- 
tor of mean values X = (ReC, ImC) used to minimize the 
Hilbert-Schmidt distance. Here we report few examples of 
the comparison between the results already obtained using (01 
with that coming from (lAU . As we will see either the two 
definitions coincide or 5' and 5 are monotone functions of 
each other. Since the definition (f2]i corresponds to an easily 
computable measure we conclude that it represents the most 
convenient choice. 

Let us first consider the Fock state g = \p) (p\. According 
to dAU . the reference Gaussian state is given by a displaced 
thermal states r' = D{C)v p D' [ (C). The overlap between g 
and t' is given by 



r c 2 

■ exp ■ 



K []p)(p\S] = -—_ 

1+p I l+pj\l+p 

c 2 



p 



X Ln 



(A2) 



The maximum of ( IA21 > is achieved for C = 0, which coincides 
with the assumptions C = Ti[a\p)(p\]. 

Let us consider the quantum state ( TTOb obtained as the so- 
lution of the loss Master Equation for an initial Fock state 
|p)(p|. The unconstrained Gaussian reference is again a dis- 
placed thermal state r' = D(C)is pv D^ (C), and the overlap is 
given by 



X Ln 



£ l+p?7 



Again, since the overlap is maximum for C = Tr[ag p (i])] = 
0, both definitions give the same results for the non- 
Gaussianity. 

Let us now consider the Schrodinger cat-like states of (0. 
The reference Gaussian state is a displaced squeezed thermal 
state, with squeezing and thermal photons as calculated be- 
fore. The optimization over the free parameter C may be done 
numerically. In Fig. |4]we show the non-Gaussianitiy, both as 
resulting from ( IA1I ) and by choosing C = Tr[ags] as in (O, 
as a function of e. The two curves are almost the same, with 
no qualitative differences. 
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FIG. 4: Non-Gaussianity of a Schrodinger cat-like state as a func- 
tion of the superposition parameter <j>, with either C obtained by nu- 
merical minimization (solid) or with C = Tr[og] (dotted). (Left): 
a = 0.5; (Right): a = 5. 



